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ABSTRACT 

In the framework of a QCD inspired relativistic potential model, we eval- 
uate radiative decay rates of heavy mesons and their coupling with soft pions. 
The agreement with the existing experimental data is satisfactory. In the limit 
itlq — > oo one obtains results in agreement with the Heavy Quark Effective 
Theory and is able to predict the values of the relevant couplings; in partic- 
ular, for the scaled B*Bir strong coupling constant g, we find that the non 
relativistic constituent quark model prediction g — 1 is modified, by the in- 
clusion of the relativistic effects due to the light quarks, to the value g = 1/3, 
in agreement with recent QCD sum rules determination. 
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In this letter we wish to analyze two problems. The first one is the study of the flavour 
conserving D* decays: 



D* -> £> 7 



D* -> Dtt . (2) 

Even though we do not posses a full set of experimental data yet, results from CLEO II 
collaboration |TJ already yield useful information on (M and (0) and provide constraints 
on the theoretical models. The second problem addressed by the present letter is the 
determination of the strong coupling go*D-K , defined by: 

< D°(k) ir + (q)\D* + (p,e) >= g D * Dn e^ (3) 

in the limit m c — ► oo. This coupling has been considered by a number of authors || |3|, |], 0. 
In fact, it is of interest in chiral effective theories of heavy mesons that describe the strong 
interactions of heavy (Qq) mesons with chiral Nambu-Goldstone bosons or light vector 
mesons, as well as their couplings to weak and electromagnetic currents H [7], R [TU], |TT ]. 



In the m c — > oo limit go'D-n can be written as follows fll"2"fl: 



2m D 

gD* Dtt = ~^g (4) 
fir 

and one can show that, in the constituent quark model g ~ 1 |TB|, |T3J] However, 
recent analyses indicate smaller values; for example from the study of the semileptonic 
D — > tt£u£ decay one estimates g = 0.4 (this value is obtained using a chiral effective 
theory and performing the limit m c — ► oo 0, Q). Also recent QCD sum rules analyses 
of the coupling (§) in the heavy quark infinite mass limit point to small values (g in the 
range 0.2 - 0.4) @. 

Prompted by these results one may wonder if this departure from the naive quark 
constituent model might be understood as a consequence of the neglect of the relativistic 



2 The value quoted in ref . [p~3| (g — 1) stems from the knowledge of the spin configurations of heavy 
mesons in the non relativistic quark model, while the authors in ref. fl4| find a slightly different value: 
g ~ 0.8, obtained in a calculation considering mock mesons (see references therein). 
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motion of the light antiquark q inside the meson D. We can test this simple explanation 
of a low value for g in a well defined model for the heavy mesons interactions, based on 
a constituent quark picture of the hadrons; the strong interaction between the quarks is 
described by a QCD inspired potential and the relativistic effects due to the kinematics 



are included in the wave equation |L5|, [16|, [17] 



Let us briefly describe the main features of this model. The heavy hadrons D a and 
D* made up by the quark Q and the antiquark q a are described by the states: 



\D a {p) >= i J dk ip(k + xp, — k + (1 — x)p) x 



v/3 

x tf(k + xp, r, a) d\(-k + (1 - x)p, s,/3)|0 > , (5) 

Sap (-e^qy, 

x b ] {k + xp,r,a) d{(-k + (1 -x)p,s,/3)\0 > , (6) 



|2?:(p,e)>= ^| ( e ^ rs J d k^k + xp,-k + (l-x)p)x 



where a and /3 are colour indices, r and s are spin indices, M an d d\ are creation operators 
of the quark Q and the antiquark q a , carrying momenta k + xp and — k + (1 — x)p 
respectively. The meson wave function ip{k + xp, —k + (1 — x)p) satisfies the Salpeter 
wave function ||18|f], which includes relativistic effects due to the kinematics explicitly: 



|y (k + xp) 2 + + y[— £ + (1 - x)p\ 2 + - y / m^~+~^|^(A; + xp, —k + (1 - x)p) 
+ / k, k')ip(k',p- k') = (7) 



which is valid in a moving frame where the meson D (or D*), having mass m^, has 
momentum p; the wave function ip is normalized as follows: 



dk\tp\ 2 = 2Jm 2 D +p 



whereas the instantaneous potential V coincides, in the meson rest frame, with the 
Richardson potential [19[] and in the r-space takes the form: 



V{r) 



8tt 



33 - 2n 



-A Ar 



/ 



/(Ar) ' 
Ar 



(9) 



1 This equation arises from the bound-state Bethe Salpeter equation by considering the instantaneous 
time approximation and restricting the Fock space to the Qq pairs; for more details see |fl6|. 
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where A is a parameter, rif is the number of flavours and []: 

n Jo q 

In order to avoid unphysical singularities |20|| , we assume that V(r), near the origin, is 
constant: 

K W = V(tm) (r < r M = ^-f ) . (11) 

The values of the parameters, as obtained by fits to meson masses, are as follows: m u = 
m d = 38 MeV; m s = 115 MeV, m c = 1452 MeV, m b = 4890 MeV, A = 397 MeV, 
A = 0.6. We note that fits of the heavy meson masses do not constrain light quark 
masses uniquely; in particular one could obtain similar results for the meson masses by 
considering , e.g. = m u ~ 100 MeV and including a small and negative constant 
term Vq in the potential. 

In order to compute the amplitudes for the decays (|l|) and (H), i.e. go'D-w in © and 
gv given by ( e=electron charge): 

< D + (k)\j;- m -\D*+(p,e) >= g v e W/3 e^V , (12) 
we have to express the currents in terms of quark operators. We write: 

1/2 

Qi(Q,r,ah^qj(^,s,/3) , (13) 



where E^q) = Jmf + q 2 , Qij is the quark charge matrix and % qj represent the usual 
quark field operators. 

We also consider the axial vector current A^, which is obtained from (|T3|) by substi- 
tuting 7 M — > 7^75 and Qij with the appropriate flavour matrix. The matrix element of 
the axial current can be written as follows: 

2 This potential grows linearly when r — > oo and follows QCD predictions for small r. 



ln(l + q 2 ) q 2 



(10) 



je.m. 



dq dq' 

W? \_E i {q)EAq')\ 
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<D°(k)\A fl \D* + (p,e)>=-i{ e M (m D . + m D )A x {q 2 



e ■ q 



mo + tud* 



(p + k)^A 2 (q 2 ) 



e-q 



2m D * gil [A 3 (q 2 ) - A (q 2 )}} 



(14) 



where 2m D *A 3 = {m D + m£ ) *)A 1 + (m D , — m£>)A 2 (q = p — k). 
Taking the derivative of A M , we obtain (J 5 = id^u): 

{m u + m d ) < D°(k)\J 5 \D*+(p, e) >= -i(e ■ q) 2m D ,A {q 2 ) . 



(15) 



For small q 2 the matrix element on the l.h.s. of (]T5p is dominated by the 7r + pole, therefore 
we obtain, for q 2 small: 



QD'Dtt 



ml 



q 2 2rriD* 



mi 



Ao(q 2 



(16) 



which gives the ttDD* coupling in the chiral limit (q 2 = 0): 



2m D * 

9d*d-it = — f — A (0) . 



(17) 



We can now compute gy in (|12|) and A (0) in ([[7]) by our model; we closely follow the 



approach described in ref. |T7| to which we refer the interested reader for further details. 
We obtain: 



where: 



9v 



Aq A q 



-i 



A 



2m 
1 

2m 



I r°° 1 / 



k' 2 



3E Q (E Q + m Q ) 
k 2 \ 
3E q (E q + m q )J 



(19) 
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where eg, e q are the heavy and light quark electric charges, Ej = Jk 2 + m 2 (j = q,Q and 



we put m u = m rf 
by the equation: 



and is related to the wave function ip of eq. (0) with p = 



u(fc) 



(20) 



£t(A;) can be obtained by solving (|7]) numerically by the Multhopp method [21]. In comput- 
ing eqs. (|18| , 19) we have put mo = vtlb,* = 2.025 GeV which is the theoretical value ob- 
tained by solving eq. (0). As for the other channels, we have put: mg = mg. = 5.33 GeV 
m D a — m D* — 2.05 GeV rriB a = ms* = 5.366 GeV. Aq and A q are the mass parameters 
whose values are reported in Table I. It may be useful to compare ( |I8D with the result of 
the constituent quark model [22], where: A q = = m u = 0.335 GeV or A q = m s = 0.45 
GeV and Aq = m c = 1.84 GeV for the charm case. 

Let us now consider gwD-n as given by eq. fllTD - Let us observe that in the present 



q 2 



[m D , 



or 



approach one obtains the value of the form factors at q 2 
(ms* —ttib) 2 ), whereas we actually need A)(0); since for small q 2 A (q 2 ) has a pole at the 
pion mass, this approximation could in principle introduce relevant differences, expecially 
for the D case. However we observe that by choosing the potential (||), we have neglected 
spin-spin terms, which corresponds to the approximation mo = m^*; for consistency we 
have therefore to put q^ = as well, which is what we shall do. Using the tecniques of 



ref. [H] we obtain the following equation for the form factor A\ in (14) 



Mlm) 



Moreover, one has: 



1 



4m d Jo 



dk\u(k)\' 



E q + m q 
E n 



k 2 



3(E q + m q ) 2 



(21) 



A 3 (q 2 



m D + m D * 



2m 



A(q 2 )' 



m D - m D > 



D* 



2m 



A 2 (q 2 



D* 



(22) 



A 3 {0) = A o (0) ; 



(23) 



therefore, in the limit mo = rap* one gets: 



A 3 (0)=A 1 (0), 
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(24) 



provided that A2(0)/mQ goes to zero when m c 
for rriQ ^> m c using relations in ref. [17 []. 

A numerical analyses gives the result: 



oo, which we have tested numerically 



A)(0) = 0.40 
A)(0) = 0.393 



(D case) 
(B case) 



(25) 
(26) 



Our results are therefore, in the chiral limit : 

9d*d-k ^ 12.3 , 
9b*Btt - 31.7 , 



(27) 
(28) 



which shows a deviation of only 2 % from the scaling result gwD^I ' 9b*b-k = mD/n^B- 



Using the results of eqs. (|T8|), (p7|), (2S) and Table I, we are able to compute the widths 
for the possible strong and radiative D* and B* decay channels. These results are reported 
in Table II together with the available experimental information. We observe that these 
results compare favourably with the experimental data whenever they are available; they 
are also in agreement with an analysis of these decays performed by us in the framework 
of HQET for mesons containing one heavy quark ||. It should be observed that in the 
our calculation we had no free parameters. 

Let us now discuss the determination of the scaled constant g defined by (HI): 



g = A (0) = A 1 (0) 



Am o Jo 



tin 



3(E q + m q ) 2 



(29) 



It is interesting to consider the non-relativistic limit, where: E q ~ m q ^> k. In this 
limit one obtains: 



1 



9 



2m D Jo 



dk\u{k)\' 



(30) 



3 We have used the analogous of eq.(36) of p7fl ; one obtains eq. ( p4| ) provided that the difference 
between xd* and xd (see (||), for their definition), does not go zero more rapidly than I/tuq. 



6 



because of the normalization condition (^). Eq. ([30]) reproduces the well known con- 
stituent quark model result |T3], pp . 

Let us now take in (^) the limit m q — > 0, which is possible since we work in the chiral 
limit and there is no restriction to the values of m q in the Salpeter equation (incidentally 
our fit of the meson masses is obtained with the rather small value m q = 38 MeV). In 
this case, we obtain: 



9 = 1 (31) 

It is worth to stress that the strong reduction of the value of g from the naive non 
relativistic quark constituent model value g = 1 (eq.(|30j))to the result ( |3~1~D has a simple 
explanation in the effect of the relativistic kinematics taken into account by the Salpeter 
equation. Furthermore, our asymptotic value of g is in the range of values (g = 0.2 — 0.4) 
determined in [[| by means of QCD sum rules in the limit of infinite heavy quark mass 
and assuming massless light quarks, that is in the same approximation that we have used 

in (grp. 

In conclusion, the potential model has allowed us to describe the radiative decays 
of heavy mesons, obtaining a rather good agreement with other similar approaches and 
with existing experimental data. Moreover, it has given us the possibility to give a simple 
explanation, based on the relativistic kinematics, for the small value of the strong coupling 
constant g of (|3]), (|j) predicted by QCD sum rules in the mg — > oo and its deviation from 
the non relativistic quark model result. 
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Table Captions 



Table I Mass parameters of eq. fllSl) (values in GeV). 



Table II Radiative and hadronic decay widths of B* and D* mesons. 
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Table I 



Decay mode 


Aq 


A g 


D* -»■ D-f 


1.57 


0.48 


D* s - Dsl 


1.58 


0.497 


B* -> £7 


4.93 


0.59 


5 S * - 5,7 


4.98 


0.66 



Table II 



Decay rate/ BR 


theory 


experiment 


T(D* + ) 




46.21 KeV 


< 131 KeV If 


BR(D* + 


-> D+7T ) 


31.3% 


30.8 ±0.4 ±0.8% 


BR{D* + 




67.7% 


68.1 ±1.0 ±1.3% 


BR(D* + 


^D+7) 


1.0% 


1.1 ± 1.4 ± 1.6% 




F(D* U ) 




41.6 KeV 




BR{D*° 


-> £)°7r°) 


50.0% 


63.6 ±2.3 ±3.3% 


BR(D*° 




50.0% 


36.4 ±2.3 ±3.3% 




T(D* S ) = 


r(L>: - ^7) 


0.382 KeF 






T(B*~) = 


= r(B*~ -> s-7) 


0.243 KeV 






T(B*°) = 


= 1(5*° -> 507) 


9.2 10~ 2 KeF 






F(B* S ) = 


r(s; - 5 s7 ) 


8.0 10~ 2 KeV 
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